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Abstract — Although memristive devices with threshold voltages 
are the norm rather than the exception in experimentally 
realizable systems, their SPICE programming is not yet common. 
Here, we show how to implement such systems in the SPICE 
environment. Specifically, we present SPICE models of a popular 
voltage-controlled memristive system specified by five different 
parameters for PSPICE and NGSPICE circuit simulators. We 
expect this implementation to find widespread use in circuits 
design and testing. 

Index Terms — Memristive devices, memristor, memristor 
model, threshold dynamics. 



I. Introduction 

IN the last few years, circuit elements with memory, 
namely, memristive (TJ, memcapacitive and meminductive 
systems have attracted considerable attention from different 
disciplines due to their capability to store and manipulate 
information on the same physical platform [3 1. However, when 
combined into complex circuits, progress in this field signifi- 
cantly relies on the available tools at our disposal. One such 
tool is the SPICE simulation environment, commonly used in 
circuit simulations and testing. While several SPICE models 
of memristive 0, 0, 0, Q, 0, (9), GUI, memcapacitive 
[6 1 and meminductive elements are already available, they 
typically [4|, [5|, |6|, |7 |, rely on physical models without 
a threshold (see, e.g., Refs. HQ, El)- 

Threshold-type switching is instead an extremely important 
common feature of memristive devices (for examples, see Ref. 
0) and, due to physical constraints, likely to be common in 
memcapacitive and meminductive elements as well. Indeed, it 
is the threshold-type switching which is responsible for non- 
volatile information storage, serves as a basis for logic oper- 
ations fl3l . fl4l . etc., and therefore, it can not be neglected. 

In the present paper we introduce a SPICE model for 
a memristive device with threshold voltage that has been 
proposed by the present authors lfT31 . Using this type of mem- 
ristive devices, we have already demonstrated and analyzed 
several electronic circuits including a learning circuit lfT31 . 
memristive neural networks [16], logic circuits [14], analog 
circuits ifTTl and circuits transforming memristive response 
into memcapacitive and meminductive ones fl8l . These pre- 
vious results thus demonstrate the range of applicability of 
the selected physical model. As a consequence, we expect its 
SPICE implementation to find numerous applications as well. 

Yu. V. Pershin is with the Department of Physics and Astronomy and USC 
Nanocenter, University of South Carolina, Columbia, SC, 29208 
e-mail: pershin@physics.sc.edu. 

M. Di Ventra is with the Department of Physics, University of California, 
San Diego, La Jolla, California 92093-0319 
e-mail: diventra@physics.ucsd.edu. 

Manuscript received May XX, 2012; revised June YY, 2012. 




Fig. 1. Sketch of the function f(V M ) for (a) a > and p > and (b) 
a = and > 0. 



II. SPICE MODEL 

The equations describing memristive systems can be formu- 
lated in the voltage- or current-controlled form 0J. I n some 
cases, a voltage-controlled memristive system can be easily 
re-formulated as a current-controlled one and vice versa 0. 
Let us then focus on voltage-controlled memristive systems 
whose general definition (for an nth-order voltage-controlled 
memristive system) is given by the following relations 

I(t) = R^(X,V M ,t)V M (t), (1) 
X = f(X,V M ,t) (2) 

where X is the vector representing n internal state variables, 
Vjv/(t) and I(t) denote the voltage and current across the 
device, and Rm is a scalar, called the memristance (for 
memory resistance). 

A specific realization of a voltage-controlled memristive 
system with threshold has been suggested by the present 
authors in Ref. lfT31l . Such a memristive system is described 




Fig. 2. Schematic of the SPICE model. The memristive device functionality 
is organized as a subcircuit consisting of a behavioral resistor R, current 
source j and capacitor C. The voltage across the capacitor (at the node x) 
defines the resistance of R. 
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Fig. 3. Memristive device directly connected to a voltage source V(t). 



by 

I = X^Vm, (3) 

^ = f(v M )[e(v M )e(R off -x) + 

9{-V M )e(X-R on )}, (4) 

with 

f(V M ) = f3V M + 0.5 (a - 0) [\V M + V t \ - \V M - V t \] , (5) 

where V t is the threshold voltage, R on and R f / are limiting 
values of the memristance Rm = X, and the ^-functions 
(step functions) are used to limit the memristance to the 
region between R on and R a ff- The important model pa- 
rameters are the coefficients a and j3 that characterize the 
rate of memristance change at \Vm\ < Vt and \Vm\ > Vt, 
respectively. These two coefficients define the slopes of the 
J(Vm) curve below and above the threshold (see Fig. [TJ. 
When a = (Fig. [TJb)), the device state changes only if 
\Vm\ > Vt. Note that Eqs. <|3j— <|5j are written in such a way 
that a positive/negative voltage applied to the top terminal with 
respect to the bottom terminal denoted by the black thick line 
always tends to increase/decrease the memristance Rm (the 
opposite convention has been used in Ref. [15|). 

The SPICE model for these devices is formulated following 
the general idea of Ref. (4). For NGSPICE circuit simulator, 
the memristive system is realized as a sub-circuit combining 
a behavioral resistor R (a resistor whose resistance can be 
specified by an expression), a current source f, and a capacitor 
C. Table [TJ presents the code of the sub-circuit. Its second line 
(Bx ...) defines the current source with the current specified 
through ternary functions. (A ternary function is defined in 
the code as a ? b : c , which means "IF a, THEN b, ELSE c" 
|[P9l .) The purpose of these functions is to limit Rm between 
R on and R ff- The third line of the code in Table [I] specifies 
the capacitor C (Cx ...) with an initial condition. The fourth 
line (Rmem ...) defines the behavioral resistor whose resistance 
takes the same numerical value as the voltage across the 
capacitor. The next line (June ...) provides the function / 
according to Eq. Q. 

For PSPICE circuit simulator, the SPICE model of mem- 
ristive device with threshold is formulated slightly differently 
without the use of behavioral resistor. Instead, we employ an 
additional current source playing the role of behavioral resistor 
[20|. In addition, in order to avoid convergence problems, 
the function / in Eq. <|5j should be smoothed. In the most 
important case of a = 0, the smoothing of / is straightforward. 
Table |II]presents the code for PSPICE circuit simulator for this 



.subekt memristor plus minus PARAMS: Ron=lK Roff=10K Rinit=5K 
alpha=0 beta=lE13 Vt=4.6 

Bx x I='(fl(V(plus)-V(minus))>0) && (V(x)<Roff) ? {fl(V(plus)- 
V(minus))}: (fl(V(plus)-V(minus))<0) && (V(x)>Ron) ? {fl(V(plus)- 
V(minus))}: {0}' 
Cx x 1 IC={Rinit} 
Rmem plus minus r={V(x)} 

.func fl(y)={beta*y+0.5*(alpha-beta)*(abs(y+Vt)-abs(y-Vt))} 
.ends 



TABLE I 

NGSPICE IMPLEMENTATION OF EQS. j5J-j5}. 



TABLE II 

PSPICE IMPLEMENTATION OF EQS. (3)-{5} FOR a = 0. 



case. In Table [II] nul and nu2 are smoothing parameters used 
in smoothed step functions /2 and /3 (although we prefer to 
use different smoothing parameters for functions of voltages 
and resistances, a common smoothing function could also be 
used). We have verified that simulation results are identical in 
both versions of SPICE. 

III. Example 

Let us consider a memristive device with threshold directly 
connected to a sinusoidal voltage source V(t) = Vq sin(27r^i) 
as presented in Fig. [3] The circuit simulations are performed 
as a transient analysis of the circuit taking into account initial 
conditions (the uic option of J ran) within the NGSPICE circuit 
simulator. In our simulations, we consider two different types 
of memristive devices with threshold corresponding to two 
cases of functions /(Vm) as presented in Fig. |TJ In the first 
case (that can be dubbed as a memristive device with a "soft" 
threshold) the coefficients a, f3 > and a < j3. In this case, 
the memristance changes at any V 7^ 0. However, the change 
is faster when the applied voltage magnitude is above the 
threshold voltage (\V\ > Vt). In the second case (Fig. [TJb)), 
a = 0. Consequently, the memristance changes only when 
the applied voltage exceeds the threshold voltage (\V\ > V t ). 
This second case is closer to the actual behavior of many 
experimentally realizable memristive systems [3|. We call this 
type of systems as memristive devices with "hard" threshold. 

Fig. [4] presents selected results of our simulations showing 
the circuit dynamics at long times (the initial transient interval 



.subekt memristor pi mn PARAMS: Ron=lK Roff=10K Rinit=5K 
beta=lE13 Vt=4.6 nul=0.0001 nu2=0.1 

Gx x value={fl(V(pl)-V(mn))*(f2(fl(V(pl)-V(mn)))*f3(Roff-V(x))+f2(- 

fl(V(pl)-V(mn)))*f3(V(x)-Ron))} 

Raux x 1E12 

Cx x 1 IC={Rinit} 

Gpm pi mn value={(V(pl)-V(mn))/V(x)} 

.tunc f 1 (y)={beta*(y-Vt)*f2(y-Vt)+beta*(y+Vt)*f2(-y-Vt)} 

.func f2(y)={l/(exp(-y/nul)+l)} 

.func f3(y)={l/(exp(-y/nu2)+l)} 

.ends 
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Fig. 4. Time evolution of the applied voltage V(t), current I(t) and memristance Rj^(t) for memristive devices with (a) soft (a = 0.1/3), and (b) hard 
(a = 0) thresholds. The vertical dashed lines - corresponding to the onset of switching - serve as guide to the eye. The simulations parameters are as follows: 
the applied voltage amplitude V = 5V, V = 0.05GHz, R on = lkfi, R oJJ = lOkfi, R M (t = 0) = 5kf2, (3 = 10 10 kn/(V s), V t = 4.6V. 



is omitted). We consider two types of memristive devices - one 
with a soft and another with hard thresholds {a = 0.1/3 and 
a = 0, respectively) - and plot the applied voltage, current 
and memristance as functions of time for these two cases at 
a frequency v — 0.05GHz. Clearly, in both cases, the current 
through the device is not of the simple sine form. The plot of 
memristance as a function of time demonstrates that the range 
of memristance change in (a) is larger than in (b) (actually, in 
(a), Rm switches between R on and R ff)- The vertical dashed 
line in Fig.|4ja) helps noticing that in Fig.[4|a) the memristance 
starts changing as soon as the sign of applied voltage changes. 
In Fig. [4jb), instead, the change of Rm occurs solely when 
| > V t . As a consequence, the shapes of Rm^) m Fig.Qa) 
and (b) are slightly different, and the steps in R]\i(t) in Fig. 
|4|b) are shifted along the horizontal axis compared to those 
in Fig. Qa). 

The current as a function of voltage at several selected 
values of v is plotted in Fig. [5] Clearly, these curves are 
typical frequency-dependent pinched hysteresis loops [TJ, 0. 
The character of the loops for memristive systems with hard 
and soft thresholds is slightly different. While for memristive 
systems with soft threshold the curve for the lowest frequency 
has the smallest loop span, the situation for the memristive 
system with hard threshold is opposite: the largest loop span 
occurs at the lowest frequency. This result, however, is not 
surprising if we take into account the fact that in the mem- 
ristive system with soft threshold the change of Rm occurs 
at lower voltages. Moreover, the insets of Fig. [3] demonstrate 



the memristance Rm (t) as a function of V(t) at a particular 
frequency. It is not difficult to notice that in the case of the 
memristive system with hard threshold (shown in the inset of 
Fig. [5jb)), Rm changes only when | V" | exceeds Vt = 4.6V . 

IV. Conclusions 

We have developed and tested a SPICE model of memristive 
devices with threshold voltage. In this model, the limiting 
conditions for the memristance are realized using ternary 
functions which adhere more closely to the actual physical 
situation, compared with the window functions approach pre- 
viously suggested [12|. The memristive device is realized as 
a sub-circuit consisting of several elements. While the present 
model is based on a single internal state variable, X, it can be 
easily generalized to more complex physical models involving 
several internal state variables. Finally, we note that threshold 
models of memcapacitive and meminductive systems can be 
implemented in the SPICE environment in a similar way. 
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